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Abstract 

In a wino LSP scenario the annihilation cross section of winos gravitationally bound in galaxies 
can be boosted by a Sommerfeld enhancement factor which arises due to the ladder of exchanged W 
bosons between the initial states. The boost factor obtained can be in the range S — 10 4 if the mass 
is close to the resonance value of M ~ 4 TeV. In this paper we show that if one takes into account 
the Sommerfeld enhancement in the relic abundance calculation then the correct relic density is 
obtained for 4 TeV wino mass due to the enhanced annihilation after their kinetic decoupling. At 
the same time the Sommerfeld enhancement in the XX ~^ W + W~ annihilation channel is sufficient 
to explain the positron flux seen in PAMELA data without significantly exceeding the observed 
antiproton signal. We also show that (e~ + e + ) and gamma ray signals are broadly compatible 
with the Fermi-LAT observations. In conclusion we show that a 4 TeV wino DM can explain the 
positron and antiproton fluxes observed by PAMELA and at the same time give a thermal relic 
abundance of CDM consistent with WMAP observations. 
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1. INTRODUCTION 



There is a great deal of activity in recent years in trying to interpret the excess of hard 
Dositron events reported by the PAMELA experiment |l| as a possible signal for dark matter 
2). In the framework of the minimal supersymmetric standard models (MSSM), one expects 
the wino to be the lightest supersymmetric particle, and hence a dark matter candidate, in 



the anoma 
models j3 : 



mediated supersymmetry breaking model 



3|-[6( as well as some string inspired 



8) . More over, such a wino dark matter is known to give the right cosmological 



relic density for a relatively heavy wino mass of 2 TeV (9]. It is also widely recognized that 
such a heavy wino dark matter can naturally account for a hard positron signal from their 
oair-annihilation into the W + W~- channel, followed by the leptonic decay of the W bosons 
10Nl3|. It is equally well known, however, that the model would need a very large boost 
factor to match the observed magnitude of the positron signal. 

One of the main theoretical mechanisms invoked to explain the above mentioned boost 



factor is the so called Sommerfeld effect 



13 



15) . It is recognized, however, that in order to 



get a large enhancement of DM annihilation cross-section (boost factor), the DM mass has 



to lie very close to one of the Sommerfeld resonances |16Hl8| , in which case the Sommerfeld 
enhancement will also have a profound effect on the DM relic density. In this work, we shall 
assume the wino DM mass to lie very close to the first Sommerfeld resonance of about 4 
TeV, and see if the Sommerfeld effect can simultaneously account for the right cosmological 
relic density as well as the large boost factor required to explain the size of the hard positron 
signal. Recently a quantitative analysis of the Sommerfeld effect in the DM relic density 
calculation has been made by Feng et al ["bol . [2^] in the context of a new physics scenario, 
where the effect arises from multiple exchanges of a light gauge boson in a hidden sector. 
We shall follow their procedure closely for incorporating the Sommerfeld effect in the relic 
density calculation for a wino DM of the MSSM, where it arises from the multiple exchanges 
of the standard W boson. We shall see below that for a wino mass of about 4 TeV, the 



Sommerfeld effect can indeed account for the right WMAP compatible relic density 21 1 
as well as a large boost factor, required to explain the PAMELA positron excess. We 
shall also see that increasing the wino mass to 4 TeV helps make the predicted anti-proton 



signal compatible with the PAMELA measurements 



22 



23] of p/p ratio by shifting the 



anti-proton peak to higher energies. The flux ratio e + / (e + + e ) measured by PAMELA 



2 



1] gets a contribution from the secondary positrons generated by cosmic ray electrons and 



protons. The measurement of the (e + + e ) flux by Fermi-LAT 24|, |25] along with the 



measurements o: 



heavier nuclei flux spectrum in the (10 — 600) GeV range by HEAO-3 



26|, ATIC-2 [23] and CREAM 



fix the parameters of primary cosmic ray background 
and the diffusion parameters in the galaxy. This leaves us essentially with the boost factor 
as the free parameter to adjust in comparing the theory with the Pamela positron data. 
Having fixed the boost factor, which in our case turns out to be S = 10 4 , we show that the 
7-rays produced by the DM annihilation do not exceed the diffuse 7-ray observations by 



Fermi-LAT 



29]. 



The work is organized as follows. In section 2, we summarize the wino DM model along 
with the procedure for incorporating the Sommerfeld effect in the DM relic density calcula- 
tion. Then in section 3, we present the results for the DM relic density after incorporating 
the Sommerfeld effect along with the corresponding boost factor as functions of the DM 
mass. In section 4 we compare the model predictions for the hard positron and antiproton 
events with the experimental data. We conclude with a summary of our results in section 5. 



2. SOMMERFELD ENHANCEMENT IN WINO DM MODEL 

In MSSM with universal gaugino mass M1/2 at the GUT scale, the one-loop renormalized 
gaugino masses at weak scale are in the ratio, 

Mi : M 2 : M 3 :: a x : a 2 : a 3 ^ 1 : 2 : 7 (1) 

and the wino is not the LSP. 

The wino LSP scenario is realized in the anomaly mediated supersymmetry breaking 



(AMSB) model , where the gaugino and scalar masses arise from supergravity break- 
ing in the hidden sector via super- Weyl anomaly contributions. The gaugino masses are 
proportional to the gravitino mass 1713/2 

= lf47r m3/2 ' = 47r m3/2 ' = 4~^ m3/2 ^ ' 

and the renormalised gaugino masses at the weak scale are in the ratio, 
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Mi : M 2 : |M 3 | :: — ai : a 2 : 3a 3 ~ 2.8 : 1 : 7.1 (3) 
5 



.3 



In the minimal AMSB model the slepton masses come out negative at the weak scale , 
and this situation is ameliorated by adding a common scalar mass term mo at tree level . 
This model has four free parameters, m 3 / 2 , tuq, tanf3 and sgn(fi). Besides the wino turns 
out to be the LSP in some string inspired models [3, S|- Both the wino annihilation and 
the Sommerfeld enhancement processes of our interest are controlled by its isospin gauge 
coupling a 2 to W boson. So our results are essentially independent of the underlying SUSY 
model; and depend only on the wino LSP mass. 

The scattering or annihilation cross section of non-relativistic particles in the initial state 
can be substantially changed if there is a long-range force between the incoming particles 
which distorts their wave function. The corrected cross section due to the distortion of 
incoming states from the plane wave can be calculated from the wave function of the two- 
body system in the attractive potential of the exchanged light particles. Taking M as the 
initial particles mass and the exchanged boson mass and a the coupling, the L=0 partial 
wave obeys the Schrodinger equation 

1 d 2 ij + a e _ m ^ = _ Mv 2^ (4) 



M dr 2 r 

with the boundary condition if)'(r) = iMvip(r) and ip(r) = e ', vr at r — > oo. The Sommerfeld 
enhancement factor is S — | C°°) 1 2 / 1 (0) 1 2 - The Sommerfeld factor can be calculated by 



solving the Schrodinger equation (j3J) numerically |l6Ml8j . In this paper we use an analytical 



approximation for the Sommerfeld factor which can be written as 



20. 



30 



3l|, 



na sinh 



12 Aft) 



S{v, M, m ) = — (5) 



r I, •<>,!,(- - i -cos(2tt ^ : ' fU/J " J 



The annihilation cross section of the S-wave initial states gets enhanced at low velocities 
by the factor S. For the W ladder processes taking = M w and a = a 2 = 1/30, the plots 
for S as a function of the DM mass M at different velocities is shown in Fig [TJ We see that 
the maximum enhancement takes place at M ~ 4TeV for any given velocity. 

For the calculation of the boost factor relevant to PAMELA signal we find with DM mass 
M = 4TeV, the peak enhancement at rms velocities 0.3 x 10~ 3 the Sommerfeld enhancement 
is S — 10 4 . It is possible that there are cold pockets in the galaxy where the rms velocities are 
lower (v ~ 10~ 4 ), and in these cold pockets the DM annihilation can get a larger boost factor 
(S ~ 10 5 ) and they can make significant contribution to the PAMELA signals [16]. In that 
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FIG. 1: Sommerfield enhancement from W exchange as a function of the DM mass for different 
relative velocities. 

case the boost factor of 10 4 represents an average DM contribution from such cold clumps 
and the smooth halo, which is parameterized here via an effective velocity of 0.3 x 10 -3 . 

For the calculation of relic density we use the relation (JSJ) for S(v, M) to calculate the 
thermal average of the cross-section which is a function of the temperature T or x = M/T, 
by averaging the cross-section over all velocities at a given temperature, 

r 3/2 roo 

(av) = —= / dv (av) S{v, M) v 2 e~ xv /4 (6) 

This thermal averaged cross section is used in the calculation of relic abundance which we 
describe in the next section. 



3. WINO RELIC DENSITY WITH SOMMERFELD FACTOR 

The thermal averaged annihilation cross section (av)(x) determines the relic density n x 
through the Boltzman equation 

dn. 



— + 3Hn x = -(av)(nl-n x l q ) (7) 

where n Xeq = g x (M 2 / (2nx)) 3 ^ 2 e~ x is the the number density of x produced by the back- 
reaction //' — > XX a t a given temperature. The Boltzman equation can be written in terms 



of dimensionless variables Y x = n x /s (where s = (2tt 2 /45)g* s T 3 is the entropy density) and 

x as, 

f = (*?<*> - *£(«>) (8) 

where 

\(x) = (^-) l/2 MM Pl (^] (av)(x) (9) 



,45/ 

where g* and g*,, are the effective degrees of freedom of the energy density and entropy 
density respectively. The freeze-out temperature is defined as the solution of 

Y(x f ) = (1 + c)Y eq (x f ) = (1 + c) (0.145) ^ x 3 / 2 e~ x f ; (10) 

9*s 

where the constant c ~ 1, g x — 2 for neutralinos; and the freeze-out temperature turns out 
to be given by Xf ~ 20. Below the freeze-out temperature Y eq (x) can be dropped from ([8]) 
and we can write the solution of OH]) in the present epoch as 

1 1 f Xs j Hx) 

+ / dx^r (11) 



Y(x s ) Y(x f ) ./,, x 2 

where x s is the temperature where the co-moving density of \ does not change noticeably 
and the integration can be terminated. The DM density fraction in the present universe is 
then given by 

n = Y (x,)soM (12) 

Pc 

where So = 2918cm -3 is the entropy of the present universe and p c = h 2 8.1 x 10 _47 GeV 4 
is the critical density. The observed value of CDM density from the seven year WMAP 
data is Q c h 2 = 0.1123 ± 0.0105 (3a) 21| (where h is the Hubble parameter in units of 100 
km/s/Mpc). 

From the calculation of relic density in the wino LSP models {9] it is known that the 
WMAP relic density Q c h 2 = 0.1123 ± 0.007(2<j) is attained if the wino mass is in the range 
M2 = (2.1 — 2.3)TeV. A wino mass of 4 TeV results in an overabundance by a factor of 3 
larger than the WMAP limit. 

In the following we calculate the relic density of WINO LSP by including the Sommerfeld 



factor in the annihilation cross section. In 



32| an approximate form of Sommerfeld factor 



was included in the annihilation cross section and it was found that wino mass which gave 
the correct relic density was in the range M 2 = (2.7 — 3.0) TeV . 
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In this paper we incorporate the Sommerfeld effect following Feng et al [19|, |20| , taking 
into account the fact that the DM particles do not share the same temperature with radiation 
bath below a the kinetic decoupling temperature but in fact cool faster than radiation. The 
Sommerfeld factor becomes very effective at temperatures below kinetic decoupling, and 
that solves the problem of DM overabundance at M=4 TeV. 



For the wino DM, the annihilation channels are x°X° ^ W + W with a cross section [32 1 



(H = -^r (is) 

and the dominant co-annihilation channel x°X~ //' with the cross section for the 5 = 



initial state 



f \ 1 71 a 2 ft A\ 

(av) = £ (14) 

2 M 2 K J 



The chargino-neutralino mass difference is AM ~ 200MeV, so the co-annilation channel will 
not take place below AM ~ (3/2)T or T = 133.3MeV. 

Kinetic decoupling: It has long been recognized that after the freezout the dark 
matter distribution at some point will not have the same temperature as the radiation bath 



but cool more rapidly 33M35I]. It has also been pointed out that in cases where the DM 
annihilation has strong velocity dependent cross sections as in the case of Sommerfeld effect 
the kinetic decoupling of DM will result in rapid annihilation of DM much below the freeze- 
out temperature 19, 20, 0, ^57\. In this section we calculate the temperature of kinetic 
decoupling of heavy wino DM and incorporate it in the wino relic density calculation. 

The xo DM will have the same temperature as the radiation bath of standard model 
particles till the rate of momentum transfer due to scattering 

T 

T k = n r (av)— (15) 

is larger than the Hubble expansion rate H = 1.66^/gTT 2 / 'Mp\. The kinetic decoupling 
temperature T^d is defined as the temperature where YkiTkd) = H(Tkd)- For T > My/ 
the scattering process XoW + — > XoW + has a cross section (av) = ^ot^jM 2 and the kinetic 
decoupling condition gives the kinetic decoupling temperature to be 

T w = 2.5 x f V2 = 5.4MeV ( —) ^ . (16) 

kd \a 2 M Pl ) V4TeVy V ; 

But at this temperature there are no in the radiation bath, so the kinetic decoupling 

temperature for this interaction channel is = (2/3)M w = 54GeV. 



Other processes that maintain temperature equality between the DM and the radiation 
bath are the elastic scattering with relativistic fermions, fx —* fx°- I n the absence of a 
Higgsino component in the DM there is no Z exchange diagram, and the elastic scattering 



34] 



process will come from s and u-channel sfermion exchanges with the cross section 

<ffl ' > = 12TO ^(MpW <17) 

where If is the isospin of / , Mj is the sfermion mass and Ef = (3/2)T is the fermion energy. 
The kinetic decoupling temperature of the relativistic fermion-wino scattering process will 
then be given by 

/(M?-M 2 ) 2 M\ 1/4 
T ^ 42 ( f M Pl ) < 18 > 

If we take M = 4TeV and Mj = lOTeV, then the kinetic decoupling for this channel takes 
place at T[ d = 5.2GeV while for M; = 5TeV,T/ d = 1.7GeV. 

Finally the quasi-elastic scattering x°f'^X ± f' w iU maintain the kinetic-coupling of x° 
with radiation down to a temperature T = (2/3) AM = 133MeV, assuming a typical charged 
and neutral wino mass difference AM = 200 MeV 0, [^l . Below the temperature 133MeV 
there will be no charginos in the heat bath so the kinetic decoupling of the winos from the 
heat bath will take place at Tkd = 133MeV. Since this is the lowest of the kinetic decoupling 
temperatures from various processes discussed above, we will take Tkd = 133MeV for the 
calculation of wino relic abundance. In Fig [2] we show the effect of the kinetic decoupling 
temperature on the relic density. Early de-couplings enhance the annihilation and lower the 
relic abundance. We take Tkd = 133MeV as the effective kinetic decoupling temperature for 
calculating the relic abundance. 

At temperatures below Tkd the temperature of the decoupled non-relativistic particles 
is related to their momenta as (3/2)T = Mv 2 /2. As the velocities of NR particles in the 
expanding universe goes down with the scale factor as v oc 1/a, DM temperature falls with 
the scale factor as 1/a 2 while the radiation temperature goes as 1/a. The integration in 
(eqn. [TTT) in the interval x > x^d = T^d/M is performed by replacing x in the integrand 
by x 2 /xkd, which corresponds to the temperature of the DM particles. This is because 
the thermal equilibrium in this region is maintained by scattering between DM particles 
via multiple W boson exchange (Sommerfeld ladder). The rapid reduction in velocities 
or temperatures of DM particles below the kinetic decoupling temperatures results in a 
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FIG. 2: Wino relic density with Sommerfeld enhancement with two different kinetic decoupling 
temperatures. T^d = 133MeV (bottom panel) is assumed for the wino relic density calculation. 
The 3<t band of WMAP relic density [21 1 is indicated by the horizontal lines. 

large enhancement in the annihilation cross section and reduction of the relic abundance in 
the late universe (T < Tkd)- With these inputs we carry out the integration in equation 
( II ip numerically to obtain Y(x s ). We plot Y(x s ) in Fig [3] and find that Y(x s ) remains 
constant after x s > 4 x 10 7 which corresponds to the temperature ~ lOOkeV. Below this 
temperature there is a chemical decoupling of the winos in the sense that their comoving 
abundance remains constant. This can be physically understood as follows. In the range 
Xkd < x < x s , the DM particles are kept in thermal equilibrium by their mutual scattering. 
But at x > x s the scattering rate falls below the Hubble expansion; so the DM particles no 
longer maintain a Maxwell-Boltzmann thermal distribution. Consequently the low velocity 



9 



particles with enhanced annihilation rate ([5} are decoupled from the high velocity particles. 
The latter escape annihilation, leading to a constant comoving abundance of DM at x > x s 
and this process is called chemical decoupling. 
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FIG. 3: Relic density as a function of cutoff temperature x s . There is a fall at x s = 10 5 due to 
kinetic decoupling, and the chemical decoupling occurs at x s = 4 x 10 . 

In Fig H] we show the relic density as a function of wino mass close to the Sommerfeld 
resonance. We see that two narrow bands of wino masses M = (3.977 — 3.983)TeV and 
M = (3.944 — 3.951)TeV are compatible with WMAP measurements at 3a. 

We choose the central values M = 3.98TeV (which is the center of the higher mass band) 
to compute the positron, electron, antiproton and diffuse 7-ray flux to compare with the 
observations of PAMELA and Fermi-LAT . 



4. COMPARISON WITH PAMELA AND FERMI-LAT DATA 

The satellite based PAMELA experiment measures the flux spectrum of e + / (e + + e~) [l| 
in the energy range of (10 — 100) GeV and p/p 22|, |23| in the range (1 — 180)GeV. It is 
found that the e + /(e + + e~) ratio exceeds the flux estimated from astrophysical sources by 
a large margin at higher energies. However the p flux is within the expected range of what 
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FIG. 4: Allowed bands of WINO mass compatible with the relic abundance measured by WMAP. 



is expected of secondary p produced from the primary CR protons from AGN's and other 
sources. Any DM annihilation or decay model must explain the large flux of positrons and 



the paucity of p's in the PAMELA signa 
and 7-rays in the by Fermi- L AT 
the Fermi-LAT data 



24 



25 



In addition there are measurements of (e + + e ) 



39 



29], and the DM signal must be consistent with 



45j . In this section we check the predictions of the heavy wino DM 



model with the PAMELA and Fermi-LAT data. 

We use the MicrOMEGAs 46[ to calculate rate of of positrons and antiprotons produced 



by the annihilation of 4TeV wino DM. As input we take the following SUSY parameters 
/I = 9 TeV, M 2 = 3.837 TeV and the other gauginos in the ratio M x : M 2 : M 3 = 2.8 : 1 : 7.1 
while the squark and slepton masses Mq = Mj = 10 TeV. We take tan/3 = 10 and // = 9TeV. 
As mentioned earlier, however, the results are insensitive to all the SUSY parameters other 
than the LSP mass. We calculate the annihilation cross section of a 3.98 TeV wino LSP 
and the branching spectrum (^) into electrons, positrons, anti-protons and photons using 
MicroOMEGAs. We fit the positron and antiproton spectrum from MicrOMEGAs with an 
analytical functions which is then input in the GALPROP code to calculate the observed 
fluxes. 
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TABLE I: Diffusion parameters and boost factor used as input in GALPROP. 



For electrons/positrons we use 

' dN P 



E7 1 



5 exp 



0.018- 



Ee 



1.5 



-0.4 



(19) 



while for protons/antiprotons the we use the analytical fit, 



dN p 
dEn 



E- 1 



5.5 



exp 



-6.35 



Ml 



0.25 



3 x 10" 



Ep 
Ml 



(20) 



and for 7-ray spectrum we use 



dN^ 
dK 



E. 



-1 



2000 



E~ 
Ml 



exp 



-5 



Ey 

Ml 



0.24 



(21) 



where E e , £/ p and £J 7 are in units of MeV, while M x r> - 



4 TeV. The integrated number 
of electrons/positrons, antiprotons and photons per wino pair annihilation are N e = 14.2 , 
N p = 1.59 and iV 7 = 26 respectively. We show the result of the output from MicrOMEGAs 
along with our analytical fits f JT9|20|) and |21] in Fig |5j 

The cross section from MicrOMEGAs of the wino annihilation to W + W~ is 5.63 x 
10 _27 cm 3 /s. We use the Galactic propagation code GALPROP 



47 



48j to obtain the 



positron and antiproton fluxes at the earth. We input a boosted cross section (<jv)qS = 
5.63 x 10~ 23 cm 3 /s in the GALPROP code along with the analytical fits of (dN/dE) (fl"9ll20l . 
We choose the isothermal dark matter distribution profile and calculate the flux of e + , e~ , p, p 
on earth arising from the wino annihilation in the galaxy. 

These inputs from MicroOmegas goes into the the source term 

of the propagation equation for the cosmic ray density (at galactic radius r and with mo- 
mentum p) 481 ]. 
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FIG. 5: Top left panel shows positron spectrum, top right panel shows antiproton spectrum and 
bottom panel shows gamma ray spectrum from wino annihilation. Light shaded curve shows the 
fitted functions (fT9ll20|) andl2Tl 

where D xx is the diffusion coefficient which describes the scattering of CR by the random 
galactic magnetic fields, V c is the convection velocity of the bulk CR in the galaxy (and 
this term represents the scattering of the CR by the background CR "wind"), D pp is the 
diffusion coefficient in momentum space (which represents acceleration in turbulent B fields), 
p is the energy loss due to radiative decay and the final two terms represent the possible 
fragmentation or radioactive decay of the CR nuclei. This equation is solved over a diffusion 
zone represented by a cylinder whose origin and axis coincide with our galactic disk. The 



half height of this cylinder is typically L ~ 1 — 10 kpc 43j. The energy dependent diffusion 
coefficient is parameterized as D xx (f,E) = Dq xx E 5 . The primary spectrum for all nuclei is 



parameterized as[43j 



ij) = -—E-^~ s (24) 
2 Dq xx 

and a similar expression for electrons with 7„ replaced by j e . For a given set of input 
parameters characterizing the primary spectra for any CR species, the GALPROP code 

13 



propagates the CR flux and gives the flux observed at the earth. We choose the set of 



parameters D 0xx , L, V c , 5, V4(Alfven velocity), ^p, 7„ (tabulated in Table]!] 



Boron/Carbon ration in CR consistent with measurements by HEAO-3 26], ATIC-2 



which gives the 



and CREAM [28j experiments as shown in FigSl This is not the only parameter set which 
gives the correct B/C spectrum, but in our scan of parameters we find that this set gives a 
smallest contribution to the p/p flux ratio. 
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FIG. 6: Boron/Carbon ratio for the assumed 



compared with data from HEAO-3 
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26(], ATIC-2 



actic diffusion parameters shown in Table 1. 



23] and CREAM 



281 ] experiments. 



The primary-electrons spectral index 5 e = 2.5 is chosen to fit the Fermi (e + + e~) data. 
The primary proton flux N p and electron flux N e are also fixed by fitting with the (e + + e~) 
and B/C observations. This set of parameters fixes the background p/p, e~ + e + and the 
7-ray photon. In Fig]7]we show the (e + + e~) signal from the 3.98 TeV wino DM model 
with observations from Fermi-LAT 2m l25 |. 

With the parameters for the background thus fixed we vary the boost factor trying to 
get a fit for the (e + /(e + + e~) data while at the same time being consistent with the p/p 
data. We find that the boost factor S = 10 4 gives a good fit to the PAMELA positron data 
with a x 2 /d.o.f = 1.03 while being consistent with the PAMELA anti-proton data (with 
X 2 /d.o.f = 1.68). 

The total flux of positrons and antiprotons from the background as above plus the DM 
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FIG. 7: The (e~ + e+) flux for the 3.98 TeV wino DM compared with FERMI-LAT data (2J, |25( 
Dashed denotes the CR background and dotted line is the DM annihilation signal. 

signal taken from GALPROP is shown in Fig [8] and Fig [9] respectively. 
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FIG. 8: Positron flux ratio for the 3.98 TeV wino DM compared with Pamela data Dashed 
line shows background from cosmic ray secondary positrons. 



The 7-ray spectrum observed at earth from DM annihilation photons is given by 

QdQ j p 2 {l)dl 

J los 



, 1 . . dAL 



(25) 
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22, 



23J. Dashed line shows background from cosmic ray secondary antiprotons while dotted line 



shows DM signal. 

where p(l) is the DM density along the line of sight (los). In addition to the prompt photons 
from DM, photons are generated also by synchrotron radiation, inverse- Compton scattering 
from the CMB and infrared starlight from the primary electrons. The prompt photons are 
sensitive to the large scale DM distribution function p(l) along the line of sight. We find that 
the Isothermal distribution gives a smaller gamma-ray flux compared to other DM profiles 
like NFW, Moore or Einasto. We plot the diffuse 7-ray flux for the M = 3.98TeV wino 
model for both the NFW , 

Pnfw = Po^ (l + — j > Po = 0.26GeV/cm 3 ,r s = 5kpc, (26) 

as well as the isothermal profiles 

l + m 2 



Piso = Po 



1+ f 



r_\2 



p Q = 0.3GeV/cm 3 , r = 8.5kpc, r s = 20kpc, 



(27) 



The local DM density for both profiles is p Q = 0.3GeV/cm 3 and r s is the size of the DM 
halo. In Fig{T0] we show the plots for the integrated 7-ray signal at large galactic latitudes 
b > 20° observed by Fermi-LAT 29j.The 7-ray signal from 3.98 TeV wino annihilation is 
consistent with measurements from Fermi-LAT. At lower energies the primary contribution 
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to the 7-ray flux is from astrophysical sources like GRB's, gamma ray pulsars and AGN's 
which is beyond the scope of the present work. 
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FIG. 10: Diffuse 7-ray flux for the 3.98 TeV wino DM compared with FERMI-LAT integrated 7-ray 
flux at galactic latitudes b > 20° [29]. Dashed line shows background. The prompt photon signal 
is smaller in the isothermal profile (top panel) compared to the NFW profile (bottom panel). At 
lower energies the primary contribution to the 7-ray flux is from astrophysical sources like GRB's, 
gamma ray pulsars and AGN's which is beyond the scope of the present work. 
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5. CONCLUSIONS 



We have investigated the possibility of reconciling the predictions of a heavy wino dark 
matter model with the WMAP data on relic density as well as the hard positron signal 
reported by the PAMELA experiment via the Sommerfeld enhancement factor. We find 
that they can be simultaneously explained if one assumes the wino mass to be very close 
the Sommerfeld resonance mass of about 4 TeV. In that case the large Sommerfeld enhance- 
ment of dark matter annihilation cross-section below the freeze-out temperature reduces the 
present dark matter relic density to the range of the WMAP measurement. Moreover the 
Sommerfeld enhancement can boost the present dark matter annihilation cross-section by a 
large factor of about 10 4 , as required to explain the size of the PAMELA positron signal. 
At the same time the large dark matter mass helps to alleviate the discrepancy with the 
PAMELA positron data by pushing up the predicted peak to still higher energies. We also 
find that the model predictions for (e + + e~) and gamma-ray signals are broadly compatible 
with the Femi-LAT observations. 
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